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Abstract 

We prove the existence of global-in-time weak solutions to a version of the parabolic-parabolic Keller-Segel system 
in one spatial dimension. If the coupling of the system is suitably weak, we prove convergence of those solutions 
to the unique equilibrium with an exponential rate. Our proofs are based on an underlying gradient flow structure 
with respect to a mixed Wasserstein-L^ distance. 

Resume 

Une note sur I’analyse variationelle du systeme de Keller-Segel parabolique-parabolique a une 
dimension spatiale. Nous prouvons I’existence de solutions faibles globales en temps d’une variante du systeme 
de Keller-Segel parabolique-parabolique a une dimension spatiale. Si le couplage du systeme est assez faible, nous 
prouvons la convergence de ces solutions vers I’equilibre univoque a une vitesse exponentielle. Nos preuves se 
reposent sur une structure de flux de gradient dans I’espace produit des espaces Wasserstein et . 


1. Introduction and main results 

We consider the following version of the Keller-Segel model for chemotaxis in one spatial dimension: 

ut{t, x) = {uxit, x) + u{t, x)Wxix) - xu{t, x)vx{t, x))x, (1) 

Vt{t, x) = Vxxit, x) - Kv(t, x) + xu{t, x), (2) 

where t > 0 and x S R, and the sought solution {u, v) is subject to the initial conditions 

u(0, x) = (x), v{0,x) = v^{x). (3) 
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We require that k > 0, x G K and the confinement potential W G C^(]R) to be bounded from below and 
to grow at most quadratically, i.e. W < W(a;) < Ax^ + B for all a; G IR and some W, A,B gR, and to 
have bounded second derivative Wxx- 

It is known that (1)&(2) possesses a variational structure since it can formally be written as a gradient 
flow of the (non-convex) entropy functional "H : X —>■ Roo (see formula (4) below) with respect to the 

compound distance dist((M, v), (u',v')) := u') + ||r' — v '\\\-2 on the space X := x 

where (^ 2 (IR),W 2 ) is the space of (absolutely continuous) probability measures - or their densities, 
respectively - on R. with finite second moment m 2 , endowed with the L^-Wasserstein distance W 2 . The 
entropy H is defined as 


TL{u, v) = 


1 /p 

I u log u + uW + -Vx + 2 ~ X^v 


+ 00 , 


dx, if / wlogudx < 00 and v G H^( 

Js. 

otherwise. 


( 4 ) 


In this note, we sketch another possible application of the method in [17] to prove the existence of weak 
solutions to (1)&(2) and to analyse their long-time behaviour. There, it has been shown that in the case 
of a porous-medium-type diffusion for u on R^, a global-in-time weak solution always exists, and that 
it converges exponentially fast to the unique equilibrium if the coupling of the system is suitably weak. 
In the one-dimensional setting at hand, some parts of the proofs simplify compared to those in [17] due 
to a gain in regularity. In contrast to that, the case of linear diffusion causes the difficulty of a missing 
time-uniform a priori estimate for u in L'"(R) for some m > 1. 

The cornerstone of our variational analysis is the so-called minimizing movement scheme (see e.g. [1,9]) 
for the construction of an approximate time-discrete gradient flow w.r.t. the distance dist: 

For each step size r > 0, let (m°,v°) := {u^,v^), and then define inductively for each n G N: 

«,v”)g argmin (^dist(('u,-u), -f'H( m, v)). (5) 

(u.ij)e.3^2(R)xL2(R) \ZT / 

Further, introduce the piecewise constant interpolation {ux ,Vr) : R+ -t ^ 2 (R) X L^(R) by 

Ur(t) = Ux, Vr{t) = Vx for alH G ((u — 1)t,ut] . (6) 

This hybrid variational principle has been exploited previously for Keller-Segel-type systems [4,5,14,16] 
in higher spatial dimensions and also in other applications, e.g. [10,12]. For the vast literature on the 
behaviour of the Keller-Segel system and its variants, we refer to the review articles by Horstmann [8] 
and Blanchet [3] and emphasize that the one-dimensional model on bounded spatial domains has been 
explicitly investigated by Osaki and Yagi [15] and Hillen and Potapov [7], leading to similar results as 
proven here. 

We obtain the following statement on the existence of global-in-time weak solutions: 

Theorem 1.1 (Existence) Assume that x> dnd W are as mentioned above and that the initial con¬ 
dition satisfies G logdx < 00 and G i7^(R). Define, for each t > 0, a discrete 

solution by means of (5)0(6). Then, there exists a vanishing sequence \ 0 (A: —>■ 00 ) such that 
{uxf,,Vxf,) converges to a weak solution {u,v) to (l)-(3) in the sense that (1) holds in the sense of distri¬ 
butions, whereas 12) and (3) hold almost everywhere. One has for all T > 0.' 

Uxf, u narrowly in the space of probability measures .^(R), pointwise with respect to t G [0,r], 

Vxk -t V in T^(R), uniformly with respect to t G [0,T], 
mg Ci/2([0,r];(^2(R),W2))nLi([0,T];L“(R))nL2([0,T];L2(M)), 

^/fiG L‘^{[0,T];H\R)), ulogu G T“([0, T]; L\R)), 
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veC°{[0,T]xR)r\H^{[0,T];L^{R))r\L°°{[0,T];H\R))r\L^{[0,T];H^{R)). 

In particular, for fixed t > 0, u{t,-) is nonnegative, continuous and bounded. The second component 
V is bounded and continuous in both variables, whereas its nonnegativity could also be obtained starting 
with a nonnegative initial condition. 

Our result on the long-time behaviour of the weak solution from Theorem 1.1 reads as follows: 
Theorem 1.2 (Convergence to equilibrium) Assume in addition to the hypotheses of Theorem 1.1 
that W is Xo-convex for some Aq > 0 and that k > 0 is strictly positive. There exist s > 0, C > 0 and 
L > 0 such that for all x = £ ^ (0, s), the following statements hold: 

(a) The system {1)&(2) possesses a unique stationary state € (t^ 2 lXL°°)(R) x II^(R) satisfying 

u°° = Ueexp{—W -|-eu°°), Ue > 0 such that ||m°°||li = 1, 


(b) One has := min(K, Aq) — eT > 0 and for all t > 0, the weak solution {u, v) to (l)-(3) from Theorem 
1.1 admits the estimate 

||M(t) - u^Wli + W2(M(t),u“) + sup \v{t) - z;“| + ||n(t) - 

a:eR 

< C(9{(u°, v°) - 9{(u°°, 

i.e. (u(t),v(t)) converges exponentially fast with rate Ag to the equilibrium as t ^ oo. 


2. Sketch of proof for Theorem 1.1 

The crucial step in the proof of Theorem 1.1 is to verify that the discrete solution (^Ut,Vt) is well- 
dehned and regular enough to allow for passage to the continuous-time limit t 0 in a strong sense. 
Once obtained, we can proceed as in [16,17] establishing an approximate weak formulation which turns 
into the weak formulation of the time-continuous equation as t 0. We prove the following 
Proposition 2.1 (Minimizing movement) For each t > 0 and {u,v) G X, the functional 

"Hr ('luj v) := ;^dist^(-, {u, v)) + % 

2t 

possesses a minimizer {u,v) G £^20^) x H^fR) with /j^uloguda; < oo. Moreover, there exist constants 
Kq, Ki, K 2 > 0 such that 

+T|kxa;||i 2 < ATo [ {ulogU - u\ogu) dx + Ki{\\v\\hi - \\v\\hi) + K 2 t(\\v\\hi + 1). (8) 

JR 

In particular, if in addition v G 1I^(R) and m log m dx < oo, then v G H'^{R), ^Ju G H^(R) and 
u G L°°(R). 

Proof. First, in one spatial dimension, there exists Co > 0 such that ||x||l“ < ||w||po4 < Co||v||//i. 
Moreover, for some Ci > 0, one has 

/ ulogwdx > —Ci(m 2 (M)-I-1)^^^. 

JR 

From this, we easily see that for all {u,v) G .^ 2 (R) x II^(R) with Jj^ulogMdx < oo, we have 

[ u\ogudx + W+ ^||x,^||i 2 - |x|Co||v||//i < 'H{u,v) < oo. 

Jr ^ 
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Using the triangle inequality for dist and Young’s inequality, we deduce coercivity of 

'Hr{u,v\u,v) > i||w||l^i + ^m 2 (u) - C. 

Thus, by the Banach-Alaoglu, Arzela-Ascoli and Prokhorov theorems, a minimizing sequence (un, Vn)neN 
for 'Hr(-lu,v) converges - at least on a subsequence - to some limit (u,v) G <^ 2 (R) x with 

/jjulogudx < oo: ^ u in Vn —>• v locally uniformly in R and Un ^ u narrowly in 

With respect to these convergences 'Ht{-\u,v) is lower semicontinuous, which is clear except for the term 
/]g UnVn dx. We employ a truncation argument similar as in [16] to prove l.s.c. for this remaining term, 
and consequently obtain the minimizing property for {u,v). It remains to prove the additional regularity 
estimate (8). We investigate the dissipation of H along the (auxiliary) 0-flow (W®,V®)s>o w.r.t. dist 
generated by the 0-geodesically convex functional 


S{u, v) 


u log u - 




dx 


on X. Elementary calculations yield, since we have V| = — kV'*: 

d 


ds 


niu%u),v%v)) < 


Using the Sobolev inequality 


-4(v^)2 + + y(V*)" 




dx. 


(9) 


we eventually arrive at 

AniW{u),V%v)) < -2\\{ViFu\l, - l||v:, - nV^Wl. + Y^^vWh + C^. (10) 

Finally, we use the flow interchange lemma [13, Thm. 3.2] to obtain S{u,v) + rD^'H(u,w) < £{u,v), 
which yields (8) in combination with (10) and lower semicontinuity as s 0. □ 

Proceeding as in [16,17], we end up with a weak solution {u, v) to (l)-(3) with the properties 

u G L“([0,r];L"(R)), u,, G L“([0, T]; e L2([0, T]; L"(R)). 


We immediately deduce that v G L°°([0,T] x R). We now show that v is continuous in both arguments. 
In fact, for all bounded intervals I CR, v belongs to the anisotropic Sobolev space 1U^’^([0,T] x I) with 

P = 1 ^ , the spectral radius of which is less than 1. Since in this case 1U^’^([0, T]xl) (g (^^([O, T] x I), 

the claim follows (for details on anisotropic spaces, see e.g. [2,11]). 


3. Sketch of proof for Theorem 1.2 

It is easily shown that the additional assumption of Ao-convexity of the confinement lU yields bounded¬ 
ness from below of the entropy TL. We thus obtain {u°°,v°°) G H L°°)(R) x i7^(R) as minimizer of TL 
similar as in [17]. Uniqueness is proved by showing strict convexity of "H as a functional on L^(]R) x L^(]R), 
which requires small coupling strength e > 0. 

Using the properties of {u°°,v°°), we observe that the entropy can be decomposed as follows into a 
convex part C (see Proposition 3.1 below) and a non-convex, but controllable part 

'H(u, v) — = L(u, v) + eLflu, v), (11) 
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where C{u, v) = Lu(u) + Cv{v), 

£u{u) := ( [ulogu — logu°° + — u°°)] dx, wiih := W — ev°°, 

JR 

£.{v)-.= \\\{v-vnAl^ + ^\\v-v-^\\l., L.{u,v):=- [ {u - u^v - dx. 

Proposition 3.1 (Properties of C) Let e be sujjiciently small. Then, the following statements hold: 

(a) There exists Mi > 0 such that the perturbed potential is Xg-convex, where Ag := Aq — Mie > 0. 

(b) The functional £„ is Xg-geodesically convex on (.^ 2 ( 1 ^)) W 2 ) and 

< £„(u) < [ u((l0gM+ Vb®)2:)^dx. 

^ Jr 

(c) The functional Cy is n-geodesically convex on L^(R) and 

f Ik - < £4^) < ^llk - - v^Wh- 

(d) There exists M 2 > 0 such that C{u, w) < (1 + M 2 s){'H{u, v) — ,v°°)). 

Actually, in one spatial dimension, the proof of part (a) simplifies dramatically compared to [17], since 

> Ao - eK|k“||Loo > Aq - eC{n{u^,v^) + 1), 

for some constant C > 0. The proof of part (d) is mainly a consequence of the Csiszdr-Kullback inequality 
(cf. [6]) Ik — u°°\Wi < CCu{u). In fact, more is true in one dimension: There exists another constant 
C" > 0 such that 

11^ - <C^[ u{{logu-\- W^)x)‘^ dx. (12) 

Jr 

The idea of proof of (12) is as follows: We distinguish the cases where the integral on the r.h.s. in (12) 
is small or large, respectively. In the former case, we can deduce from that an L°° bound on u leading 
to the desired result using the Taylor expansion of the integrand in £„ at Uodx). The latter case can be 
treated by a suitable Sobolev interpolation in one spatial dimension. 

We now prove the central estimate leading to Theorem 1.2: 

Proposition 3.2 (Exponential estimate for £) Let {uf,vf)nGN be a family of time-discrete approx¬ 
imations obtained by (5) which converges to a weak solution (u,v) as r \ 0. Then, there exist e > 0 and 
L > 0 such that for all e G (0, e) and n G N, one has 

£«,<) < il+M2s)in{u°,v°)-n{u^,v°°))il + 2A,T)-^, (13) 

with Ag := min(Ao, k) — Ls > 0. 

Once proven, this result yields exponential convergence of L{u(t),v{t)) to zero for t —>■ 00 after passage 
to the continuous-time limit r 0. From this, Theorem 1.2 clearly follows. 

Proof. We investigate the dissipation of TL along the (auxiliary) min(Ag, «:)-flow (W^,V®)s>o of the 
min(Ag, K)-geodesically convex functional £ on X, which is associated to the evolution system 

u: = VI = (V* - n“),, - «(V^ - n“). 

First, by elementary calculations, we obtain using the decomposition (11) 

^n{W{u), V\v)) < (I - 1) + VFkx)" dx + I - v^)l dx 

+ |liz^* -+ (I -1) ll(v* --^“)lli- 
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The third term can be controlled by the first one using (12), whereas the second term is to be controlled 
by the fourth term using the inequality ||? 7 /x^||ii < C'||? 7 || which is valid in one spatial dimension. 

Taking into account the properties of C from Proposition 3.1, we end up with 

--^'H{W{u),V\v)) > 2(1 - eM) min(Ae, k)L{W{u), V%v)), (14) 

ds 

for some constant M > 0 if e is sufficiently small. The application of the flow interchange lemma [13, Thm. 
3.2] eventually yields with (14): [1 + 2t(1 — erM) min(A£, k)]£(w", v") < £(w"“^, v"“^). By iteration of 
this estimate and Proposition 3.1(d), the desired estimate (13) follows. □ 
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